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TRANSONIC  SHOCK  -  TURBULENT  BOUNDARY  LAYER 
INTERACTION  AND  INCIPIENT  SEPARATION  ON  CURVED  SURFACES 


George  R.  Inger 
University  of  Colorado 
Boulder,  Colorado 


Abstract 


A  detailed  analysis  is  made  of  weak  normal 
shock  -  turbulent  boundary  layer  interactions  on 
longitudinally-curved  surfaces  for  the  case  of  non¬ 
separating  steady  2-D  flow.  It  is  shown  that  the 
interactive  viscous  displacement  effect  on  the  lo¬ 
cal  outer  inviscid  transonic  flow  eliminates  the 
well-known  singularity  pertaining  to  a  curved  wall. 
The  inner  interaction  solution  within  the  boundary 
layer  reveals  that  curvature  moderately  influences 
the  interaction  through  the  turbulent  eddy  viscos¬ 
ity.  A  non-asymptotic  triple-deck  solution  valid 
over  a  wide  range  of  practical  Reynolds  numbers  is 
given  which  incorporates  this  effect,  and  example 
numerical  results  are  presented  and  verified  by 
comparison  with  experimental  data.  Small  amounts 
of  curvature  .01- 0&  are  found  to  moder¬ 

ately  spread  out  and  thicken  the  interaction  zone 
while  also  delaying  slightly  the  onset  of  any  in¬ 
cipient  separation  that  occurs  under  the  shock. 


Nomenclature 

Speed  of  sound  at  sonic  conditions 

Skin  friction  coefficient  (=  2t  /o  U 

“  *1  el 

Pressure  coefficient  (=  2p'/p  U  ) 

el  el 

Functions  defined  in  Eqs.  20,23 
Body  shape  function,  see  Fig.  7 
Incompressible  form  factor  (=  /  0^*) 

Wall  curvature  (=  1/R  ) 

D 

Distance  to  shock  location 
Mach  number 
Static  pressure 

Interactive  pressure  perturbation  (=  p  -  p^) 
Pressure  jump  across  incident  shock 
Wall  radius  of  curvature  (~  K 
Reynolds  number  based  on  length  £ 

Reynolds  number  (*  p  U  L  /  u  ) 
el  el  el 

Transformed  dependent  variables,  Eq.  17 
Absolute  temperature 

Velocity  components  in  x,y  directions,  re¬ 
spective  ly 

Streamwise  and  normal  interactive  disturb¬ 
ance  velocity  components,  respectively 
Undisturbed  incoming  boundary  layer  veloc¬ 
ity  in  x-direction 
Resultant  flow  velocity 

Coordinates  along  and  normal  to  body  surface 
Stretched  independent  variables,  Eqs.  18,21 
Effective  wall  shift  seen  by  interactive  in¬ 
viscid  flow 
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Shock  wave  angle.  Fig.  6 

Specific  heat  ratio 

Boundary  layer  thickness 

Boundary  layer  displacement  thickness 

Kinematic  turbulent  eddy  viscosity 

Non-dimensional  interactive  displacement 
thickness 

Resultant  flow  direction  angle 
Momentum  thickness 
Ordinary  coefficient  of  viscosity 
p/p 

Density 

Similarity  parameter  (see  Eq.  18) 

Total  shear  stress 

Basic  interactive  wall-turbulence  parameter 
Viscosity-temperature  dependence  exponent. 


Sub-  and  Superscripts 

0  Undisturbed  incoming  boundary  layer  proper¬ 

ties 

1  Inviscid  flow  conditions  ahead  of  shock 

2  Inviscid  flow  conditions  behind  shock 

e  Conditions  at  boundary  layer  edge 

inv  Inviscid  disturbance  solution  value 

w  Conditions  at  wall  surface 

(  )*  Denotes  perturbation  from  "0"  state 

1.  Introduction 

An  understanding  of  transonic  shock  -  turbulent 
boundary  layer  interactions  is  important  in  the 
aerodynamic  design  of  high-speed  aircraft  wings 
(both  ordinary  and  circulation-controlled) ,  turbine 
and  cascade  blades  in  turbomachinery,  and  air  breath¬ 
ing  engine  inlets  and  diffusors.  Since  these  ap¬ 
plications  often  involve  curved  surfaces  and  since 
a  singularity  is  associated  with  a  normal  shock  on 
a  curved  surface  in  purely  inviscid  flow,  the  in¬ 
fluence  of  wall  curvature  on  transonic  shock  - 
boundary  layer  interaction  is  an  important  basic 
and  practical  question.  Oswatitsch  and  Zierep^ 
studied  the  related  problem  of  normal  shock  im¬ 
pingement  on  a  curved  wall  in  a  purely  invisc id 
flow;  they  found  that  convex  curvature  introduces 
a  logarithmic  singularity  in  the  wall  pressure  in 
the  form  of  a  sharp  post-shock  expansion  (Fig.  1), 
and  this  phenomenon  has  indeed  been  observed  in  in¬ 
viscid  numerical  solutions  by  Emmons 2 ,  Jameson3  and 
others.  However,  in  real  flows  even  at  high  Rey¬ 
nolds  numbers,  it  is  known  that  viscous  smearing 
and  upstream  influence  effects  arise  due  to  the 
thin  boundary  layer  along  the  surface  which 
profoundly  influence  the  physics  of  the  local  in¬ 
teractive  flow'4»,).  Consequently,  in  view  of  the 
well-known  role  often  played  by  viscosity  in  elim¬ 
inating  singularities,  one  would  expect  the  bound¬ 
ary  layer  to  significantly  alter  and  perhaps  even 
eliminate  this  inviscid  curvature-induced  singu¬ 
larity.  In  the  present  paper,  we  address  this 
issue  for  the  case  of  a  weak  nearly-normal  transon¬ 
ic  shock  interacting  with  a  2-D  non-separating  tur¬ 
bulent  boundary  layer.  Our  approach  is  to  extend 
the  successful  transonic  shork  -  boundary  layer  in- 
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teraction  theory  of  Inger6  to  Include  the  new  fea¬ 
tures  introduced  by  the  combined  Inviscid  and  tur¬ 
bulent  boundary  layer  flows  on  a  curved  wall. 

Following  brief  overviews  of  the  general  physi¬ 
cal  features  of  the  problem  and  the  non-asymptot ic 
triple-deck  solution  method  in  the  flat  surface 
case,  we  first  made  a  detailed  analysis  of  the  lo¬ 
cal  mixed  transonic  inviscid  flow  structure  outside 
the  boundary  layer  when  it  is  interactively-coupled 
with  the  attendant  thickness  response  of  the  layer. 
The  results  show  that  the  inviscid  singularity  dis¬ 
appears  (the  external  flow  being  now  quite  regular) 
when  the  boundary  layer  displacement  effect  is 
taken  into  account.  We  then  examine  the  inner  dis¬ 
turbance  problem  within  the  boundary  layer,  where 
it  is  shown  that  surface  curvature  does  exert  an 
influence  on  the  interaction,  but  from  a  hitherto- 
overlooked  source:  its  effect  on  the  turbulent 
eddy  viscosity.  Using  an  approximate  incorpora¬ 
tion  of  this  effect  into  the  triple-deck  solution, 
example  numerical  results  are  given  and  validated 
by  experimental  comparisons.  A  re-interpretation 
of  curved  wall  data  in  the  light  of  the  present 
findings  is  also  given  along  with  parametric  study 
results  establishing  the  conditions  for  incipient 
separation  on  curved  as  well  as  flat  surfaces. 


2.  Typical  Features  of  the  Interaction 

As  a  background,  consider  the  typical  principal 
physical  features  of  a  non-separating  normal  shock 
-  turbulent  boundary  layer  interaction.  To  fix 
ideas,  we  consider  an  undisturbed  turbulent  bound¬ 
ary  layer  in  the  practical  Reynolds  number  range 
105  <  ReL  <  103  that  is  perturbed  by  a  normal  or 
nearly-normal  shock  weak  enough  to  avoid  separa¬ 
tion  (Mi  £1.3);  whether  the  wall  is  flat  or  curved, 
it  is  an  experimental  fact  under  these  conditions 
that  the  resulting  interaction  flow  pattern  is  of 
the  relatively  simple  type  illustrated  in  Fig.  2A. 
The  boundary  layer  spreads  out  the  shock  perturba¬ 
tion  upstream  as  well  as  downstream  to  give  a  wall 
pressure  distribution  such  as  shown  in  Fig.  2H, 
with  a  corresponding  significant  growth  of  the 
boundary  layer  displacement  thickness  particularly 
behind  the  shock,  so  that  even  though  the  original 
undisturbed  boundary  layer  is  thin  and  negligibly 
affects  the  overlying  inviscid  flow,  this  is  no 
longer  the  case  in  the  interaction  zone  astride  the 
shock  whether  or  not  the  wall  is  flat  (Figs.  2B-E). 

Thus,  the  local  outer  inviscid  flow  "sees"  the 
overall  viscous  effects  in  the  interaction  as  an 
effective  thickening  of  the  geometrical  surface  by 
an  amount  equal  to  the  interactive  displacement 
thickness  growth.  In  the  vicinity  of  the  shock, 
this  is  seen  to  involve  (a)  a  rapid  growth  in  the 
effective  body  slope  at  the  shock  foot  followed  by 
(b)  an  overall  downstream  increase  of  the  effective 
body  thickness.  At  the  very  least,  these  two  fea¬ 
tures  must  be  included  in  any  realistic  mathemati¬ 
cal  model  of  the  true  outer  inviscid  flow  in  such 
interactions.  Indeed,  th^s  conclusion  is  strongly 
supported  by  recent  studies  of  interaction  effects 
by  YoBhihara18  and  Murman18. 

3.  Non-Asymptot ic  Triple-Deck  Treatment 
of  Flat-Surface  Interactions 

3.1)  Triple-Deck  Structure 

It  is  well-known  experimentally  that  when  sepa¬ 


ration  occurs,  the  disturbance  flow  pattern  associ¬ 
ated  with  normal  shock  -  boundary  layer  interaction 
is  a  very  complicated  one  involving  a  bifurcated 
shock  pattern7,  whereas  the  unseparated  case  per¬ 
taining  to  turbulent  boundary  layers  up  to  Mjf  1.3 
has  instead  the  much  simpler  type  of  interaction 
pattern  (Fig.  2A)  which  is  more  amenable  to  ana¬ 
lytical  treatment.  With  some  judicious  simplifica¬ 
tions,  it  is  possible  to  construct  a  fundamentally- 
based  approximate  theory  of  this  flow  problem  for 
flat  walls.  For  purposes  of  orientation  and  com¬ 
pleteness,  a  brief  outline  of  the  essential  fea¬ 
tures  of  this  theory  will  now  be  given  (see  Ref.  6 
for  full  details). 

Consider  a  known  adiabatic  boundary  layer  pro¬ 
file  Mq  (y)  subjected  to  small  transonic  disturb¬ 
ances  due  to  an  impinging  weak  and  nearly  normal 
shock.  At  high  Reynolds  numbers  it  has  been  estab¬ 
lished8*10  that  the  resulting  interactive  perturba¬ 
tion  field  in  the  neighborhood  of  the  shock  organ¬ 
izes  itself  into  three  basic  layers  or  "decks" 

(Fig.  3):  (1)  an  outer  region  of  potential  invis¬ 

cid  flow  above  the  boundary  layer,  containing  the 
incident  shock  and  interactive  wave  systems;  (2)  an 
intermediate  deck  of  frozen  shear  stress-rotational 
inviscid  disturbance  flow  occupying  the  outer  90Z 
or  more  of  the  incoming  boundary  layer  thickness; 

(3)  an  inner  shear-disturbance  sublayer  adjacent  to 
the  wall  which  accounts  for  the  interactive  skin 
friction  perturbations  (and  hence  any  possible  in¬ 
cipient  separation)  plus  most  of  the  upstream  in¬ 
fluence  of  the  interaction.  The  "forcing  function" 
of  the  problem  here  is  thus  impressed  by  the  outer 
deck  upon  the  boundary  layer;  the  middle  deck  cou¬ 
ples  this  to  the  response  of  the  inner  deck  but  in 
so  doing  can  itself  modify  the  disturbance  field, 
while  the  slow  viscous  flow  in  the  thin  inner  deck 
reacts  strongly  to  the  pressure  gradient  disturb¬ 
ances  imposed  by  these  overlying  decks.  In  the 
practical  Reynolds  number  range  of  interest  here, 
we  analyze  this  disturbance  field  with  a  non- 
asymptotic  triple-deck  treatment  patterned  in  many 
ways  after  Lighthill's  approach4  because  of  its  es¬ 
sential  soundness  and  adaptability  to  further  im¬ 
provement,  its  similarity  to  related  multiple-deck 
approaches  that  have  proven  highly  successful  in 
treating  turbulent  boundary  layer  response  to 
strong  known  adverse  pressure  gradients,  and  the 
large  body  of  data  that  supports  the  predicted  re¬ 
sults  in  a  variety  of  specific  problems. 

3.2)  Formulation  of  the  Disturbance  Problem 

Outer  Potent i a 1  FI ow  Reg i on 

If  the  incident  shock  and  its  reflection  system 
are  weak  with  isentropic  non-hypersonic  flow,  we 
have  a  small  disturbance  potential  inviscid  motion 
in  the  undisturbed  uniform  flow  U  ,  M  : 


where  the  third  term  within  the  square  brackets  is 

significant  in  the  transonic  case  I'M  <1.05-1.10 

o 

e 

and  includes  shock  jump  conditions  to  this  order  of 
approximation.  Since  various  solution  methods  are 
available  in  either  transonic  or  purely  supersonic 
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flow  (in  which  case  Eq.  1  reduces  to  an  Ackeret 
problem),  we  assume  that  such  a  solution  may  be 
carried  out  for  all  x  on  y  >  60  subject  to  the 
usual  far-field  conditions  as  y  -*■  <».  The  remain¬ 
ing  disturbance  boundary  condition  along  y  =  6q 

then  couples  this  solution  to  the  underlying 
double-deck:  it  requires  that  both  v*  / U  and  p1 

be  continuous  there.  ° 

Middle  Rotational-Disturbance  Flow  Deck 


This  layer  contributes  to  and  transmits  the 
displacement  effect,  contains  the  boundary  layer 
lateral  pressure  gradient  due  to  the  interaction 
and  carries  the  influence  of  the  incoming  boundary 
layer  profile  shape.  Our  analysis  rests  on  the  key 
simplifying  assumption  that  for  non-separating  in¬ 
teractions  the  turbulent  Reynolds  shear  stress 
changes  are  small  enough  to  have  a  negligible  back 
effect  on  the  mean  flow  properties  along  the  inter¬ 
action  zone;  hence,  they  are  "frozen”  along  each 
streamline  at  the  appropriate  value  in  the  undis¬ 
turbed  incoming  boundary  layer.  This  approximation 
is  well-supported  by  detailed  experimental  stud- 
les 1 1 » 12.  Thus  the  disturbance  field  is  one  of 
small  rotational  inviscid  perturbation  of  the  in¬ 
coming  non-uniform  turbulent  boundary  layer  profile 
governed  by 
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where  UQ(y),  P  <y)  are  arbitrary  functions  of  y 

with  i  ,  &  and  t  as  constants.  Now  Eq.  3  is  a 

o  o  w 

slight  generalization  of  Lighthill's  well-known 
pressure  perturbation  equation  for  non-uniform 
flows,  which  includes  a  non-linear  correction  term 
for  transonic  effects  including  the  diffracted  im¬ 
pinging  shock  above  the  sonic  level  of  the  incoming 
boundary  layer  profile.  Excluding  the  hypersonic 
regime,  Eqs.  2  and  3  therefore  apply  to  a  wide 
range  of  initially  supersonic  external  flow  condi¬ 
tions  and  across  the  boundary  layer  except  at  the 
singular  point  Mq  *  0  (which  we  avoid  by  considera¬ 
tion  of  the  inner  deck  as  shown  below).  Whatever 
the  method  used  to  solve  this  middle  deck  disturb¬ 
ance  problem,  we  imagine  that  it  provides  the  dis¬ 
turbance  pressure  distribution  p'(x,y);  then  y- 
integration  of  Eq.  2  gives  the  disturbance  stream¬ 
line  slope  as 
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where  y 


>0  is  the  effective  wall  height  of 


eff 

the  inner  deck  defined  such  that  the  inviscid 


v'(x,y  )  and  hence  3p'/<*y(x,y  )  both  vanish 

“eff  “eff 

(see  below). 


The  corresponding  displacement  thickness  growth 
along  the  interaction  is  then  given  by  streamwise- 
quadrature  of  the  continuity  equation  integral  as 


The  Inner  Shear- Disturbance  Layejr 

This  thin  inner  deck  contains  the  significant 
viscous  and  turbulent  shear  stress  disturbances  due 
to  the  interaction.  It  lies  well  within  the  Law  of 
the  Wall  region  of  the  incoming  boundary  layer  and 
below  the  sonic  level  of  the  profile.  The  original 
work  of  Lighthill4  and  others  further  neglected  the 
turbulent  stresses  and  considered  only  the  laminar 
sublayer  effect;  while  this  yields  an  elegant  ana¬ 
lytical  solution,  the  results  are  in  error  at  high 
Reynolds  numbers  and  cannot  explain  the  ultimate 
asymptotic  behavior  pertaining  to  the  Re  00  lim¬ 
it.  The  present  theory  remedies  this  by  Including 
the  entire  Law  of  the  Wall  region  turbulent  stress- 
effects.  Note  that  our  consideration  of  the  en¬ 
tire  Law  of  the  Wall  combined  with  the  use  of  ef¬ 
fective  inviscid  wall  concept  to  the  inner  deck 
displacement  effect  eliminates  the  "blending  lay¬ 
er"3  that  is  otherwise  required  to  match  the  dis¬ 
turbance  field  in  the  laminar  sublayer  region  with 
the  middle  inviscid  deck,  since  we  impose  a  boundary 
condition  of  vanishing  total  shear  disturbance  at 
the  outer  edge  of  the  deck. 

To  facilitate  a  tractable  theory,  we  retain 
only  the  main  physical  effects  by  introducing  the 
following  simplifying  assumptions:  (1)  the  incom¬ 
ing  boundary  layer  Law  of  the  Wall  region  is  char¬ 
acterized  by  a  constant  total  (laminar  plus  turbu¬ 
lent  eddy)  shear  stress  and  a  Van  Driest-Cebeci 
type  of  damped  eddy  viscosity  model.  (2)  For  weak 
incident  shocks,  the  sublayer  disturbances  are 
small  perturbations  upon  the  incoming  boundary  lay¬ 
er;  however,  all  the  physically-important  effects 
of  streamwise*  pressure  gradient,  streamwise  and 
vertical  acceleration,  and  both  laminar  and  turbu¬ 
lent  disturbance  stresses  are  retained.  Moreover, 
since  the  form  of  the  resulting  set  of  linear  equa¬ 
tions  is  in  fact  unaltered  by  non-linear  effects, 
the  the  quantitative  accuracy  is  expected  to  be 
good  until  close  to  separation.  (3)  For  adiabatic 
flows  at  low-to-moderate  external  Mach  numbers,  the 
undisturbed  and  perturbation  flow  Mach  numbers  are 
both  quite  small  within  the  shear  disturbance  sub¬ 
layer;  consequently,  the  influence  of  the  density 
perturbations  may  be  neglected,  while  the  corres¬ 
ponding  compressibility  effect  on  the  undisturbed 
profile  is  treated  by  using  incompressible  rela¬ 
tions  based  on  wall  recovery  temperature  properties. 
(A)  The  turbulent  fluctuations  and  the  small  inter¬ 
action  disturbances  are  uncorrelated.  (5)  The 
thinness  of  the  inner  deck  allows  neglect  of  its 

lateral  pressure  gradient:  p*  ■  p'(x). 

w 

Under  these  assumptions,  the  disturbance  field 
is  governed  by  the  continuity  and  momentum  equa- 
t  ions 
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where  is  the  kinematic  eddy  viscosity  perturba¬ 

tion.  The  corresponding  undisturbed  turbulent 
boundary  layer  Law  of  the  Wall  profile  U0(y) 
governed  by 
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where  the  Van  Driest-Cebeci  eddy  viscosity  model 
/ — - 1 

with  y*  =  /t  /  p  /  u  yields  for  non-separ- 
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ating  flow  that 
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where  we  take  the  commonly-accepted  value 


(9a) 

(9b) 


A  =  26, 


We  seek  to  solve  Eqs.  6-9  subject  to  the  imper¬ 
meable  wall  no-slip  conditions  Uo(0)  =  u'(x,o) 

=  v' (x,o)  «  0  plus  an  initial  condition  requiring 
all  interactive  disturbances  to  vanish  far  upstream. 
Furthermore,  sufficiently  far  from  the  wall,  u' 
must  pass  over  to  the  inviscid  solution  uj  v  along 
the  bottom  of  the  middle  deck  governed  by  Eq,  6  plus 
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while  the  corresponding  total  shear  disturbance 
(  3u'  /  3y)  vanishes  to  a  desired  accuracy. 


Now  differentiate  Eq.  7  w. r.t.  x,  substitute 
Eq.  6  so  as  to  eliminate  u'  and  then  differentiate 
the  result  w.r.t.  y  so  as  to  eliminate  p^  ;  this 
yields  the  following  fourth-order  equation  for  v': 
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along  with  3  v1  /  3y  -  0  (i.e.,  vanishing  total 
disturbance  shear).  Once  this  v'(x,y)  field  is 
obtained,  the  attendant  streamwise  velocity  and 
hence  the  disturbance  shear  stress)  may  then  be 
found  from  Eq.  6. 

An  important  feature  here  is  the  "effective  in¬ 
viscid  wall"  position  (or  displacement  thickness) 
that  emerges  from  the  asymptotic  behavior  of  v'-far 
from  the  wall4  (Fig.  4A),  This  is  defined  by  the 
value  Yweff  where  the  "back  projection"  of  the  v^nv 

solution  vanishes;  physically,  yu  thus  represents 

e  1 1 

the  total  mass  defect  height  due  to  the  shear  stress 
perturbation  field  and  hence  the  effective  wall  po¬ 
sition  seen  by  the  overlying  inviscid  middle  deck 
disturbance  flow.  As  indicated  in  Fig.  4B,  this 
couples  the  inner-  and  middle-deck  solutions  by  pro¬ 
viding  the  non-singular  inner  equivalent  slip-flow 
boundary  conditions 

3P '  '  3y  (yeff>  =  vinv  <y«eff)  =  0  at  Uo  (ywef f }  >0 
for  the  middle-deck  solution. 


3.3)  Approximate  Solution  by  Operational  Methods 


An  approximate  analytic  solution  is  further  a- 
chieved  by  assuming  small  linearized  disturbances 
ahead  of  and  behind  the  nonlinear  shock  jump*  plus 
an  approximate  treatment  of  the  detailed  shock 
structure  within  the  boundary  layer,  which  give  ac¬ 
curate  predictions  for  all  the  properties  of  engin¬ 
eering  interest  when  M  >  1.05.  The  resulting  equa¬ 
tions  can  be  solved  by  operational  methods6*" 
yielding  the  interactive  pressure  rise,  displace¬ 
ment  thickness  growth  and  the  skin  friction  behavior 
downstream  and  upstream  of  the  shock  foot.  In  par¬ 
ticular,  the  matching  of  the  outer  two  decks  with 
the  inner  shear-disturbance  deck  in  connection  with 
the  Fourier  inversion  process  yields  the  determina¬ 
tion  of  the  upstream  influence  distance  eu,  the 
inner  deck  displacement  thickness 


H(T) 


(14) 

p-l/3 


This  equation  contains  a  three-fold  influence  of 
the  turbulent  flow:  the  profile  UQ(y),  its  curva¬ 
ture  d2U0  /  dy2  (non-zero  outside  the  laminar  sub¬ 
layer)  and  a  new  eddy  disturbance  stress  term  2  . 

o 

Eq.  11  is  to  be  solved  together  with  Eqs.  8,9 
and  the  wall  boundary  conditions  v* (x,o)  * 

3v*  /  3y(x,o)  *  0.  A  third  condition  involving  v* 
is  obtained  by  satisfying  the  x-momentum  equation 
right  at  the  wall: 

,3  ,  ,  d2p' 

~  (x.o)  -  -  (2  V  r1  (12) 

~  J  o.Z 

dy  w  dx 

The  fourth  boundary  condition  is  the  outer  inviscid 
matching  from  Eq.  10: 


and  the  interactive  skin  friction  relation 
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This  nonlinear  shock  jump  provision  plus  the  vari¬ 
ous  non-uniform  viscous  flow  effects  within  the 
boundary  layer  reduces  the  lower  Mach  number  limit 
otherwise  pertaining  to  the  linearized  supersonic 
theory  in  purely  Inviscid  potential  uniform  flow. 


(16b) 
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where  Km^n  =  ^  i-s  8*-ven  in  Ref.  6  while  the 

functions  H(T)  and  S(T),  shown  here  in  Fig.  5, 
represent  the  wall  turbulence  effect  on  the  inter¬ 
active  displacement  effect  and  skin  friction,  re¬ 
spectively.  Fig.  5  is  a  central  result,  providing 
a  unified  account  of  the  entire  Reynolds  number 
range  in  terms  of  the  single  new  turbulent  interac¬ 
tion  parameter  T  from  the  limiting  behavior  of 
negligible  wall  turbulence  effect  pertaining  to  the 
T  -►  0  limit4* 9  to  the  opposite  extreme  of  wall  tur¬ 
bulence-dominated  behavior  at  T  >>  1  pertaining  to 
an  asymptotic-type  of  theory®**0  at  very  large  Rey¬ 
nolds  numbers  where  the  inner  deck  thickness  and 
its  disturbance  field  become  vanishingly  small. 

A  computer  program  has  been  constructed  to  car¬ 
ry  out  the  foregoing  solution  method;  it  involves 
the  middle-deck  disturbance  pressure  solution  cou¬ 
pled  to  the  inner  deck  by  means  of  the  effective 
wall  shift  (Eq.  14)  combined  with  an  upstream  in¬ 
fluence  solution  subroutine.  The  corresponding  lo¬ 
cal  total  interactive  displacement  thickness  growth 
and  skin  friction  are  obtained  from  Eqs.  5  and  15, 
respectively.  If  desired,  the  attendant  boundary 
layer  shape  factor  change  along  the  interaction  may 
then  also  be  calculated  as  H =  [6*0+A6* (x) 3/0* (x) 
with  0*  given  by  an  x-wise  integration  of  the  over¬ 
all  momentum  integral  equation  for  the  total  local 
boundary  layer  since  p(x),  6*  and  Cf  are  known.  The 
incoming  turbulent  boundary  layer  is  treated  by  the 
compressible  version  of  a  universal  composite  Law 
of  the  Wall  -  Law  of  the  Wake  model  due  to  Walz13 
that  not  only  has  a  convenient  analytical  form  (see 
Appendix)  but  also  provides  a  very  general  fundamen¬ 
tal  description  of  this  boundary  layer  in  terms  of 
three  arbitrary  parameters:  preshock  Mach  number, 
boundary  layer  displacement  thickness  Reynolds  num¬ 
ber,  and  the  incompressible  shape  factor  .  This 

enables  us  to  treat  the  important  but  heretofore- 
neglected  influence  of  the  upstream  flow  history 
(pressure  gradient,  suction,  etc.)  on  the  interac¬ 
tion;  indeed,  surface  curvature  effects  on  the  pre¬ 
interactive  flow  may  thus  be  taken  into  account  al¬ 
so. 


4.  Extension  to  Curved  Walls 

Our  analysis  of  the  curved  surface  case  is  a 
composite  treatment  involving  a  basic  hodograph 
analysis  of  the  outer  inviscid  transonic  flow  due 
to  Sobieczky17  interactively  linked  to  a  wall  cur¬ 
vature  -  corrected  extension  of  the  foregoing  non- 
asymptotic  triple-deck  theory  for  the  inner  bound¬ 
ary  layer  region. 

4.1)  The  Outer  Inviscid  Disturbance  Flow 

Consider  the  local  mixed  inviscid  transonic 
flow  in  the  neighborhood  of  the  shock  outside  the 
boundary  layer,  assuming  the  shock  weak  enough  to 
be  isentropic.  Then  if  W  and  are  the  magnitude 
and  direction  angle,  respectively,  of  the  resultant 
flow  velocity  and  one  introduces  the  transformed 
dependent  variables 


s-±i  <«'>m 

(17a) 

t  *  a-1  0 

(17b) 

as  functions  of  the  stretched  independent 

variables 

X  *■  x 

(18a) 
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Y  =  j  (y+l)o  y 

(18b) 

where  the  constant  a  is  chosen  later  as  an  appropri¬ 
ate  similarity  parameter,  then  the  governing  tran¬ 
sonic  flow  equations  are: 

~l/3  as  3t 

3X  '  3Y 

(19a) 

c-1/3  3S  -  St  „ 
b  )Y  !X  " 

(19b) 

This  pair  of  equations  in  general  has  either  super¬ 
sonic  wave-type  solutions  (S  >  0,  minus  sign  in  19b) 
or  subsonic  Cauchy  -  Riemann  -  type  solutions  (S  <  Q 
positive  sign).  If  we  further  restrict  attention 
to  small  perturbations  about  the  shock  condition, 
Eqs.  19  become  linear  in  the  leading  approximation. 

Now  choose  o  such  that  upstream  of  the  shock 
S  =  =  1;  then  linearization  of  Eq.  19  gives  a 

wave  equation  with  solutions  of  the  form*7 

s  =  si  +  F1(f.)  +  CL(n)  (20a) 


The  aforementioned  solution  method  captures  all 
the  essential  global  features  of  the  mixed  transon¬ 
ic  viscous  interaction  flow,  including  lateral 
pressure  gradient  effects  and  interactive  skin  fric¬ 
tion  up  to  incipient  separation  for  a  very  arbi¬ 
trary  input  turbulent  boundary  layer  profile.  This 
has  been  verified  by  numerous  detailed  comparisons 
with  experiment  over  a  wide  range  of  Mach  -  Reynolds 
number  -  shape  factor  conditions*4* Thus  the 
resulting  predictions  (such  as  typified  by  Fig.  2) 
are  believed  to  give  a  sound  engineering  account  of 
non-separating  interaction  regions  in  practical 
transonic  flow  fields  and  provide  an  appropriate 
framework  for  extension  to  the  case  of  moderately- 
curved  surfaces. 


t  =  tx  +  Fx(0  -  Cl  (n)  (20b) 

where  £  =  X  +  Y  and  n  -  X  -  Y.  Linearizing  the 
subsonic  equations  for  perturbations  of  the  post¬ 
shock  condition  S2  <  0,  t2  yields  the  following 

solution  of  the  Cauchy-Riemann  equations17: 

Z  =  X  +  iS21/3Y  =  S2  +  it2  F2(Z)  (21) 

A  particular  case  of  a  pair  of  local  solutions 
(Eqs.  20  and  21)  pertains  to  the  flow  ahead  of  and 
behind  an  oblique  shock  of  the  "strong  branch" 
type,  i.e.,  with  subsonic  flow  downstream.  The 


transonic  shock  jump  conditions  link  the  pre-  and 
post-shock  conditions  as  follows: 


At  =  t2-  V  2'3/2  •  3  •  (1  -  S22/3)1/2(l  +  S22/ J)  (22a) 
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Thus,  if  at  some  point  P  there  is  a  supersonic  flow 
W[  >  a*  at  inclination  i?[  that  is  decelerated  by  an 
oblique  shock  to  some  W2  <  a*  at  angle  t?2  (Fig-  6A), 
then  W) ,  ^ 1 ,  W2,  $2  are  related  in  terms  of  Sj,  tj, 
S2>  t2  by  means  of  Eqs.  17a,  17b,  22a.  This  is 
the  transonic  shock  polar  (Fig.  6B),  and  the  wave 
angle  C  corresponding  to  the  flow  deflection  angle* 
is  found  from  Eq.  22b.  Then  with  given  shock  jump 
conditions  at  P,  the  local  solutions  (Eqs.  20  and 
21)  may  be  used  to  determine  the  velocity  gradients 
as  well  as  the  curvature  of  the  shock  shape  in  the 
vicinity  of  P. 


This  local  transonic  inviscid  perturbation 
field  now  can  be  linked  to  the  inner  viscous  dis¬ 
turbance  flow  by  noting  from  above  that  the  outer 
flow  sees  an  effective  body  consisting  of  the  bare 
wall  hw  "  hi  •  plus  a  slope  jump  C2  at  the 

shock  foot  followed  by  a  further  interactive  dis¬ 
placement  thickness  growth,  giving  h2  “  C„X 

f  D2X2  (see  Fig.  7).  Upstream  of  the  shocf,  an 'Ar¬ 
bitrary  linear  variation  of  inviscid  flow  Mach  num¬ 
ber  along  the  wall  is  assumed  to  allow  for  either 
decelerating,  accelerating  or  uniform  incoming 
supersonic  flow.  Now  in  terms  of  this  flow  model, 
Sobieczky17  carried  out  a  local  solution  by  means 
of  a  double  series  expansion  analysis**  of  the  up¬ 
stream  and  downstream  flow  properties  about  the 
shock-foot  origin  using  the  following  particular  so¬ 
lutions  of  Eqs.  20  and  21: 


F1  -  V 

(23a) 

G1  =  V 

(23b) 

=  A2Z  +  B2Z  In  Z 

(23c) 

where  the  downstream  (complex)  coefficients  A2,  B2 
are  related  to  their  upstream  (real)  counterparts 
and  the  shock  shape  parameters  by  the  transonic 
shock  polar.  The  results  of  this  analysis  yield 
two  basically  different  types  of  physical  behavior, 
one  regular  and  the  other  singular,  depending  on 
whether  the  interactive  deflection  effect  is  in¬ 
cluded  (C2  >  0,  D2  4  C^)  or  neglected  a  priori 

(C2  »  0,  D2  =  0^),  as  follows: 

c2  >  0,  d2  4  cx 

REGULAR  SOLUTION  INCLUDING  6*-EFFECT 
Parabolic  Oblique  Shock  Shape:  Xgh » Qy + Py2  (24a) 


where  this  solution  Implies  that  the  displacement 
thickness  effect  causes  the  inviscid  shock  to  be¬ 
come  curved  near  the  boundary  layer  edge  with  a 
slightly  oblique  angle  at  the  foot,  as  indeed  ob¬ 
served  experimentally 'S  1  e'2  1  (see,  e.g..  Fig.  8). 


SINGULAR  SOLUTION  NEGLECTING  6*-EFFECT 
Curved  Shock  Normal  to  Wall:  Xsj,  =  y2 (Q  +  R  In  y ) 

(25a) 

Logarithmically  Singular  Shocked  Flow: 

M2  =  -K2  (1  -  E2  X -F2X  lnX)  (25b) 

In  Eqs.  24  and  25  the  parameters  Q,  P,  K2,  L2>  R, 

K|,  E2,  F2  are  constant  coefficients  given  in 
terms  of  Cj,  C2  and  D2  (see  Ref.  17  for  details). 
Now  the  "bare"  curved  wall  solution  (Eqs.  25)  is  in 
fact  the  Oswatitsch-Zierep  singular  solution1  men¬ 
tioned  earlier.  It  has  been  of  use  mainly  to  veri¬ 
fy  computational  results  for  inviscid  flows  with 
recompression  shocks:  fully-conservative  codes  for 
airfoils  usually  give  results  with  sharp  post-shock 
expansion  on  the  curved  suction  surface,  demonstrat¬ 
ing  the  locally-exact  treatment  of  the  recompres¬ 
sion  shock  if  the  computational  grid  is  fine  enough1 
—  see  Fig.  9.  On  the  other  hand,  when  the  presence 
of  the  turbulent  boundary  layer  is  taken  into  ac¬ 
count  with  its  resulting  ramp-like  interactive  dis¬ 
placement  effect,  the  regular  solution  (Eqs.  24)  is 
clearly  the  appropriate  inviscid  model  even  on  a 
curved  wall  because  it  takes  into  account  the  flow 
deflection  and  oblique  shock  distortion  due  to  this 
interaction-induced  boundary  layer  thickening.  This 
major  conclusion  has  been  forcefully  corroborated 
by  the  supercritical  airfoil  flow  field  studies 
carried  out  by  Murman  et  al19:  as  clearly  illus¬ 
trated  in  Fig.  10,  when  they  include  the  influence 
of  the  interaction  "bump"  (modeled  by  a  ramp  simi¬ 
lar  to  the  one  described  here)  on  the  effective  in¬ 
viscid  solution  wall  shape,  the  post-shock  singu¬ 
larity  otherwise  occurring  on  the  curved  surface 
completely  disappears  from  their  inviscid  flow 
field  calculations. 

The  foregoing  discussion  shows,  at  least  for 
non-separating  flows,  that  the  interactive  dis¬ 
placement  thickness  effect  does  indeed  eliminate 
the  inviscid  singularity,  giving  in  fact  a  well-be¬ 
haved  external  inviscid  disturbance  flow.  Moreover, 
it  is  seen  that  this  regular  influence  of  curvature 
on  the  inviscid  flow  via  the  displacement  effect  is 
necessarily  of  order  S*/Rb  and  hence  numerically 
very  small  in  most  practical  applications.  To  be 
sure,  wall  curvature  does  have  some  influence,  as 
shown  below,  but  this  is  also  quite  regular  and 
does  not  derive  significantly  from  the  inviscid 
part  of  the  flow. 


Linear  Downstream  Acceleration:  M2-I  •  K2+L2X  (24b) 


It  is  assumed  that  the  deflection  angle  -  dj 

lies  below  the  maximum  value  of  the  shock  polar, 
which  is  reasonable  in  practice  if  the  Incident 
shock  is  weak  enough  to  avoid  separation. 

** 

The  realistic  assumption  is  implied  that  the  ra¬ 
dius  of  any  wall  curvature  is  large  compared  with 
a  typical  boundary  layer  displacement  thickness. 


4.2)  Inner  Boundary  Layer  Region 

Now  wall  curvature  influences  the  inner  solu¬ 
tion  in  two  ways:  (1)  it  alters  the  incoming  un¬ 
disturbed  turbulent  boundary  layer  profile  upon 
•  which  the  perturbation  solution  depends;  (2)  it 
introduces  new  explicit  terms  in  the  small  dis¬ 
turbance  equations.  The  first  is  by  far  the  more 
important  and  is  discussed  further  below;  the  sec¬ 
ondary  effects  (2)  are  in  fact  negligible  for  the 
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small  relative  curvatures  (6  /R  <  .02)  encoun- 

O  D 

tered  in  practice  according  to  the  following  con¬ 
siderations.  As  regards  the  (regular)  influence  of 
curvature  on  the  inviscid  outer  boundary  conditions 
via  the  displacement  effect,  the  foregoing  analysis 
shows  that  this  is  of  order  6  */R  and  hence  very 

O  D 

small  in  the  leading  approximation.  Likewise,  un¬ 
der  the  continued  assumptions  that  the  interactive 
disturbances  remain  uncorrelated  with  the  back¬ 
ground  turbulent  fluctuation  field  and  that  the  vis¬ 
cous  disturbance  sublayer  lies  within  the  Law  of 
the  Wall  region  of  the  turbulent  boundary  layer, 
the  new  terms  in  both  the  rotational  inviscid  and 
viscous  disturbance  sublayer  equations  can  be 
shown  to  be  of  the  same  order  or  smaller  and  hence 
also  negligibly  small-  Thus,  to  a  consistent  de¬ 
gree  of  first  approximation,  the  form  of  the  equa¬ 
tions  governing  the  disturbance  flow  in  the  bound¬ 
ary  layer  is  not  altered  by  wall  curvature. 

Regarding  the  basic  turbulent  flow  profile, 
once  again  the  explicit  new  curvature  terms  in  the 
flow  equations  governing  this  profile  (including 
the  centrifugal  lateral  pressure  gradient  effect) 
have  all  been  shown  to  be  of  the  order  and 

hence  negligible  (e.g.,  see  Ref.  23  and  more  recent 
work24).  However,  the  influence  on  the  eddy  vis¬ 
cosity  relation  is  known  to  have  an  order  of  magni¬ 
tude  larger  effect  (10  to  20  times  6  /R„)  on  the 

o  B 

overall  skin  friction  and  shape  factor  Hf;  since  it 
has  been  firmly  established  that  the  role  of  the 
profile  in  the  interaction  mainly  derives  from 
these  properties16*25,  this  influence  is  therefore 
deemed  the  cause  of  any  significant  curvature  ef¬ 
fect  on  the  interaction  solution. 


By  virtue  of  the  foregoing  arguments,  the  cur¬ 
vature  effect  can  be  reasonably  estimated  by  incor¬ 
porating  into  the  incoming  boundary  layer  profile 
model  appropriate  corrections  for  the  influence  of 
curvature  on  CfQ,  the  shape  factor  and  thickness. 

Examination  of  Bradshaw  s  comprehens  ive  study2  3  shows 
that  over  a  wide  range  of  parametric  and  local  con¬ 
ditions,  these  corrections  for  small  longitudinal 
curvature  and  non-separating  flows  can  be  adequate¬ 
ly  represented  by  the  engineering  approximations: 


(26b) 


where  to  this  first  order  accuracy  the  correspond¬ 
ing  effect  on  60  is  much  smaller  and  therefore  neg¬ 
lected.  Note  that  the  typical  value  6  /R_  =  .01 

o  B 

yields  a  reduction  and  increase  in  Cf  and  .  of 

1  o  ll 

10%  and  5%,  respectively,  which  are  an  order  of 
magnitude  larger  than  the  explicit  new  curvature 
terms  in  either  the  mean  or  perturbation  flow  e- 
quations  (including  the  outer  inviscid  disturb¬ 
ances).  The  use  of  Eqs.  26  with  6  /R_  as  a  fourth 

o  B 

input  parameter  in  the  aforementioned  non-asymp- 
totic  triple-deck  theory  thus  enables  a  straight¬ 
forward  appraisal  of  the  first  order  curvature  ef¬ 
fects  on  the  non-separating  normal  shock  -  turbu¬ 
lent  boundary  layer  interaction  problem.  In  addi¬ 
tion,  this  theory  has  been  generalized21*22  to  in¬ 


clude  the  influence  of  the  shock  obliquity  (provid¬ 
ed  the  post-shock  state  remains  subsonic)  which  at¬ 
tends  the  displacment  effect  on  the  external  invis¬ 
cid  flow. 


5.  Discussion  of  Results 
5.1)  Influence  of  Wall  Curvature 

A  systematic  study  of  the  wall  curvature  effect 
based  on  the  foregoing  extended  interaction  theory 
has  been  made  over  a  wide  range  of  conditions  for 
K<5  =  6  /R  ratios  up  to  the  largest  values  (  *.02) 

O  O  D  — 

normally  found  in  practical  applications.  Repre¬ 
sentative  results,  which  are  representative  of  all 
the  cases  examined  ( 1 0 3  <  Re6*  106,  1.1  <^M  <^1.30, 
1.30  £  llj,  1.7)  are  presented  in  Figs.  11  through 
13  and  will  now  be  discussed. 

The  typical  interactive  pressure  distributions 
presented  in  Figs.  11A  and  11B  show  that  wall  cur¬ 
vature  slightly  spreads  out  the  interaction,  weak¬ 
ening  the  adverse  pressure  gradient,  and  that  this 
effect  derives  primarily  from  the  increased  shape 
factor  (Fig.  11A):  The  results  for  the  pressure 
disturbance  along  the  outer  edge  of  the  boundary 
layer  (Fig.  11B)  further  bring  out  several  inter¬ 
esting  fundamental  points.  First,  a  small  subson¬ 
ic  expansion  region  is  predicted  to  occur  right  be¬ 
hind  the  shock  regardless  of  the  wall  curvature; 
this  is  an  inherent  feature  of  the  mixed  transonic 
viscous  interaction  flow  along  a  flat  surface9  as 
further  confirmed  by  detailed  numerical  solutions8 
(Fig.  12)  and  Gadd's  pipe  flow  experiments  (Fig. 
13).  Thus,  contrary  to  what  is  sometimes  alleged, 
wall  curvature  per  se  is  not  the  cause  of  such  ex¬ 
pansion  regions  in  non-separating  turbulent  inter¬ 
actions.  Second,  when  streamwise  distance  is  prop¬ 
erly  rescaled  so  that  x  -  0  (*Q)  as  is  done  in  Fig. 
11B ,  these  regions  are  perfectly  regular  as  pre¬ 
dicted  by  Eqs.  24.  Third,  although  not  the  cause, 
convex  wall  curvature  does  strengthen  the  expansion 
and  reduce  the  preceding  local  shock  pressure  jump. 
We  re-emphasize,  however,  that  this  is  a  result  of 
the  curvature  effect  on  the  boundary  layer  eddy 
viscosity  and  has  nothing  whatsoever  to  do  with  any 
inviscid  curvature  singularity. 

Figs.  14A-D  illustrate  as  a  group  the  modest 
curvature  effects  on  all  the  main  streamwise  prop¬ 
erties  along  an  interaction  for  a  set  of  condi¬ 
tions  quite  typical  of  a  full  scale  wing.  Since 
curvature  increases  and  hence  reduces  the  incom¬ 
ing  profile  "fullness,"  it  acts  to  spread  out  the 
wall  pressure  disturbance,  increase  the  upstream 
influence,  reduce  the  downstream  pressure  level  and 
(Fig.  14C)  thicken  the  downstream  boundary  layer. 
The  corresponding  influence  on  skin  friction  (Fig. 
14D)  is  also  noteworthy:  although  curvature  re¬ 
duces  the  upstream  level,  it  slightly  increases  the 
local  Cf  near  the  shock  foot  owing  to  the  reduced 
adverse  interactive  pressure  gradient  and  hence  de¬ 
lays  the  onset  of  separation.  This  prediction  is 
concordant  with  Gadd’s  assertion  (p.  32  of  Ref.  26) 
that  "there  is  some  direct  evidence  that  curvature 
affects  the  way  in  which  boundary  layer  velocity 
profiles  respond  to  the  rise  of  pressure...  ;  ...  by 
convex  surface  curvature  .  .  .  thus  separation 
might  be  delayed." 


3.2)  Comparisons  with  Experiment 


Ackeret,  Feldmann  and  Rott's  famed  experimen- 
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tal  study7  of  shock  -  boundary  layer  interaction  on 
a  plate  and  wall  in  the  choked  transonic  flow  of  a 
slightly-curved  wind  tunnel  nozzle  provides  some 
examples  of  unseparated  turbulent  flow  suitable  for 
at  least  approximate  comparison  with  the  present 
theory.  We  have  chosen  those  for  which  both  wall 
and  inviscid  pressure  distributions,  as  well  as 
displacement  thickness,  are  given.  It  should  be 
noted,  however,  that  direct  comparisons  with  these 
data  involve  numerous  uncertainties  in  converting 
to  the  theoretical  variables  (or  vice-versa):  (a) 
the  boundary  layer  thickness  and  hence  the  inviscid 
flow  edge  is  only  approximately  defined;  (b)  the 
shock  wave  location  and  shape  are  uncertain  to  with¬ 
in  .25  50  to  .50  (c)  reading  the  curves  where 

they  change  rapidly  introduces  inherent  error;  (d) 
a  significant  background  inviscid  pressure  gradient 
is  present  which  beclouds  interpretation  of  the 
cuter  fringes  of  the  interaction  zone  and  the  char¬ 
acter  of  the  incoming  "undisturbed"  turbulent  layer 
profile;  (e)  the  upstream  boundary  layer  history  is 
only  partially  understood,  especially  following 
forced  transition  cases,  and  this  together  with  (d) 
cannot  be  fully  accounted  for  in  the  local  interac¬ 
tion  model;  (f)  some  amount  of  channel  blockage  ef¬ 
fect  occurs  from  the  interactive  boundary  layer 
thickening,  which  reduces  the  post-shock  channel 
area  and  hence  the  effective  theoretical  shock 
strength  and  downstream  interaction  pressure  level 
(this  was  independently  identified  by  both  Panaras 
and  Inger25  and  Mel  nick  and  Grossman27  who  devised 
different  but  equivalent  correction  methods  for  it). 
For  these  reasons  it  is  understood  that  the  follow¬ 
ing  comparisons  are  primarily  of  qualitative  value. 

A  typical  non-separating  interactive  pressure 
field  measured  by  AF  &  R  is  illustrated  in  Fig.  15A; 
experimental  pressure  distributions  along  both  the 
surface  and  the  approximate  boundary  layer  edge  are 
compared  in  Fig.  15B  with  our  theoretical  predic¬ 
tion  (corrected  for  the  estimated  interactive 
blockage  effect*  using  the  Panaras-Inger25  method), 
while  the  corresponding  displacement  thickness  and 
skin  friction  distributions  are  shown  in  Figs.  15C 
and  15D,  respectively.  We  note  that  the  relevant 
non-dimensional  curvature  parameter  for  these  ex¬ 
periments  was  rather  small  (K  <5Q  -  .0063)  so  that 
the  predicted  curvature  effect  on  the  interaction 
is  only  slight,  as  indicated.  Regarding  the  skin 
friction  comparison  shown  in  Fig.  15D,  we  note  that 
the  "experimental"  values  were  actually  inferred 
from  measured  velocity  profiles  along  the  interac¬ 
tion  (hence  0*  and  •$*)  by  means  of  the  2-D  momen¬ 
tum  integral  equation;  considering  the  uncertain¬ 
ties  of  the  experimental  set  up  and  this  method, 
combined  with  the  present  theory's  own  limita¬ 
tions  ,  the  agreement  is  considered  good  as  regards 
both  the  magnitude  and  shape  of  the  Cf  curve.  In¬ 
deed,  in  view  of  the  aforementioned  difficulties 
in  interpreting  the  data  plus  the  approximations  of 
the  present  small  disturbance  theory,  the  overall 
agreement  is  deemed  to  be  quite  good.  In  particu¬ 
lar,  the  following  definitive  features  of  the  in- 


*  * 

For  the  significant  6  -  increases  (100-200%) 

typically  encountered  even  in  non-separating  inter¬ 
actions,  this  blockage  effect  is  found  to  be  sig¬ 
nificant  when  all  the  wall  surfaces  are  taken  into 
account;  in  particular,  for  the  relatively  narrow 
channel  dimensions  of  the  AF  &  R  tests,  this  meth- 
od  indicates  that  ^ef f  =  ' 8”  x  ^ankine-Hugoniot 


teraction  predicted  by  the  the  theory  are  well-cor¬ 
roborated:  (1)  the  magnitude,  sign  and  streamwise 

extent  of  the  lateral  pressure  gradient  effect  both 
ahead  of  and  behind  the  shock;  (2)  the  existence  of 
a  long  slow  interactive  pressure  rise  (algebraic 
rather  than  exponential)  downstream  of  the  shock; 
(3)  the  overall  streamwise  scale  and  upstream  in¬ 
fluence  distance;  (4)  the  magnitude  and  shape  of 
the  interactive  displacement  thickness  growth; 

(5)  the  local  inviscid  pressure  jump  across  the 
shock  at  the  boundary  layer  edge;  (6)  a  non-singu¬ 
lar  inviscid  subsonic  expansion  region  behind  the 
shock  due  to  the  viscous-inviscid  interaction  and 
not  surface  curvature  (note  that  the  zero  curvature 
in  this  region  is  actually  closer  to  the  data), 
consistent  with  the  regular  inviscid  post-shock 
model . 


A  second  set  of  comparisons,  illustrated  in 
Fig.  16,  involves  some  recent  DFVLR-AVA(Go)  data 
obtained  on  two  supercritical  airfoils  with  non¬ 
separating  interaction  zones  on  the  curved  upper 
surfaces.  Here,  the  theoretical  prediction  method 
used  incorporates  the  present  analysis  as  a  local 
interactive  module  within  a  global  viscous-inviscid 
prediction  program  for  supercr it ical  airfoils22.  It 
is  seen  that  the  theory  yields  predictions  of  both 
the  wall  pressure  and  displacement  thickness  dis¬ 
tributions  along  the  interaction  region  that  are  in 
excellent  agreement  with  the  measured  values  (in¬ 
deed  this  was  found  true  in  many  other  cases  not 
shown  here28).  Good  agreement  is  also  shown  in  the 
corresponding  local  skin  friction  behavior,  where 
we  note  that  the  experimental  Cf  values  were  in¬ 
ferred  from  various  streamwise  station  boundary  lay¬ 
er  profile  surveys  by  means  of  the  Ludwig  Tillman 
relationship 


Cf  =  .246(T*/T er796  RV-268  e-1-561  Hi 


(27) 


where  T*/T  =  1  +  .14  M  2  for  y  =1.40  and  P  =  .7 
e  e  r 

on  an  adiabatic  wall. 


5.3)  Incipient  Separation 


The  particular  attention  paid  by  the  present 
theory  to  the  accurate  analysis  of  local  interac¬ 
tive  skin  friction  behavior,  including  the  role  of 
wall  curvature,  makes  it  possible  to  establish  con¬ 
ditions  under  which  incipient  separation  (C.  ■*  0) 

local 


occurs . 


Thus  setting  t  -  1'  +  i 

w  wo  i, 


=  0.  Kq.  12 


yields  the  following  explicit  analytical  criterion 
for  the  onset  of  separation: 

12/3 


C’ 

P, 


2  r  3/2/< 

p  min 


f 


cc'  )3/2dx 
p«. 


/Cf  / 6  *  2/3 


(28) 


where  it  is  re-emphasized  that  C'  here  is  the  lo- 

Pw 

cal  interactive  distribution.  Kq.  28  bears  a  gen¬ 
eral  resemblance  to  a  Stratford-type29  of  incipient 
separation  relation  for  turbulent  flow,  except  that 
the  present  formula  contains  the  integrated  history 
effect  along  the  interaction  whereas  Stratford's 
result  involves  purely  local  properties  of  Cp  and 
dCp^/dx.  It  is  understood,  of  course,  that  the 
present  theory  actually  breaks  down  approaching 
such  separation  owing  to  the  combination  of  its 


8 


linearization  assumptions  and  the  Van  Dr iest /Cebeci 
wall  turbulence  model  used;  nevertheless,  Eq.  28 
does  give  at  least  an  approximate  indication  of 
where  this  will  occur  and  indeed  does  so  without 
containing  any  adjustable  empirical  constants.  Note 
that  according  to  the  present  theory  wall  curvature 
effects  on  incipient  separation  are  implicitly  ac¬ 
counted  for  via  their  influence  on  the  values  of 
CfQ  and  used  in  evaluating  Eq.  28. 

Based  on  the  foregoing,  a  systematic  parametric 
study  of  the  critical  shock  strength  for  incipient 
separation  vs.  Reynolds  number,  shape  factor  and 
degree  of  wall  curvature  was  carried  out  to  estab¬ 
lish  a  fundamental  "incipient  separation  curve"; 
the  results  are  shown  in  Fig.  17.  Also  indicated 
is  the  approximate  experimental  boundary  deter¬ 
mined  by  a  careful  examination  of  all  available 
transonic  interaction  data3^,  plus  Nussdorfer * s3 1 
well-known  Mf a? 1 . 30  criterion  for  turbulent  flow. 

It  is  seen  that  the  theoretical  prediction  of  a 

gradual  increase  in  the  Mi .  .  value  with 

1 incip.sep. 

Reynolds  number  is  in  agreement  with  the  experimen¬ 
tal  trend;  moreover,  the  theoretical  prediction  of 
only  a  small  influence  of  shape  factor  on  the  in¬ 
cipient  separation  conditions  is  also  borne  out  by 
the  lack  of  any  consistent  H-effect  for  the  same  Re 
discernible  in  the  data.  Stariewsky7®  and  Squire32 
have  experimentally  observed  similar  insensitivi¬ 
ties  to  in  transonic  and  purely  supersonic  flow 
interactions,  respectively.  The  absolute  values  of 
M^incip  sep  predicted  by  the  present  interaction 

theory  are  seen  to  be  consistently  slightly  lower 
than  the  average  experimental  value;  this  is  at¬ 
tributable  to  the  combined  effects  of  the  linear¬ 
ized  inner  deck  theory  (which  overpredicts  the 
pressure  gradient  effect  on  Cf  and  hence  too  small 
an  incipient  separation  shock  strength)  and  the  as¬ 
sumption  of  a  normal  shock  when  in  fact  most  experi¬ 
ments  likely  entail  some  shock  obliquity  (which  al¬ 
so  delays  separation  to  somewhat  higher  shock- 
strengths).  In  conjunction  with  these  results  it 
is  interesting  to  note  that  Nussdorfer' s  original 
incipient  separation  criterion,  based  as  it  was  on 
a  very  limited  base,  does  roughly  go  through  the 
average  of  the  data  although  it  does  not  account 
for  the  proper  Reynolds  number  effect. 

As  shown  in  Fig.  1 7B,  wall  curvature  in  the 
range  0  *•  •  .02  has  only  a  small  effect  on  in¬ 

cipient  separation:  at  Reynolds  numbers  below  a- 
bout  ReL  ~  5  x  107,  it  delays  the  onset  of  separa¬ 
tion  to  slightly  higher  shock  strength  for  a  given 
Hj owing  to  the  reduced  interactive  pressure  gra¬ 
dient  that  lowers  Cf^  (see  Fir.  7),  whereas  at 
high  Reynolds  numbers  Ehe  effect  ’everses  sign 
owing  to  the  eventual  dominating  influence  on  the 
curvature  solution  of  the  reduction  in  Cf  with  in¬ 
creasing  Re^.  Qualitatively,  therefore,  ?he  influ¬ 
ence  of  curvature  is  similar  to  (but  even  smaller 
than)  that  of  increasing  shape  factor  (see  Fig.  17A). 
Indeed,  from  a  practical  standpoint,  the  predicted 
curvature  effect  clearly  lies  within  the  indicated 
uncertainty  band  of  the  data  establishing  incipient 
separat ion . 

6.  Concluding  Remarks 

This  investigation  has  sought  to  delineate  the 
essential  influence  of  surface  curvature  on  non¬ 
separating  turbulent  boundary  layer  -  transonic 
shock  interactions,  to  place  in  true  perspective 
the  inviscid  curvature  singularity,  and  to  devise 


L 


an  approximate  analysis  including  curvature  that 
enables  parametric  study,  comparisons  with  experi¬ 
ment  and  an  approximate  prediction  of  incipient 
separation.  The  resulting  theory  is  an  approximate 
non-asymptot ic  triple  deck  model  of  the  interaction 
pertaining  to  the  Reynolds  number  range  10rj<Re^<108 

The  major  results  of  the  study  are  the  follow¬ 
ing.  (1)  When  the  viscous  displacement  effect  on 
the  inviscid  flow  is  accounted  for  (as  it  should 
by  very  definition  of  an  interaction  problem),  the 
surface  curvature  effect  on  non-separating  interac¬ 
tions  is  non-singular,  moderate  and  derives  from 
its  influence  on  the  eddy  viscosity  within  the 
boundary  layer.  (2)  The  region  of  sharp  pressure 
rise  followed  by  expansion  observed  along  the  edge 
of  the  boundary  layer  may  be  physically  interpreted 
as  inherent  features  of  the  viscous  mixed  transonic 
interaction  that  occur  even  on  a  flat  surface  and 
which  are  only  secondarily-influenced  by  wall  cur¬ 
vature.  (3)  Owing  to  the  slight  spreading  out  and 
hence  weakening  of  the  adverse  interactive  pressure 
gradient  it  causes,  wall  curvature  delays 
the  onset  of  incipient  separation  to  slightly  high¬ 
er  shock  Mach  numbers;  for  ordinary  practical  cases 
(K6C  <  .02),  however,  this  effect  may  be  neglected. 

Three  areas  of  future  application  warranting 
further  study  appear  to  be  of  practical  interest. 

(a)  Extension  of  the  present  interaction  theory  to 
the  unsteady  case  (examining  first  the  validity  of 
the  quasi-steady  approximation)  in  order  to  study 
unsteady  air  loads  due  to  flutter  at  transonic 
speeds.  (b)  Adaptation  to  three  dimensional  flow 
fields  of  finite-span  wings,  at  least  outside  wing / 
fuselage  juncture  or  tip  -  influence  regions.  (c) 
More  detailed  study  of  the  effects  of  shock  -  bound¬ 
ary  layer  interactions  on  transonic  internal  flows 
within  engine  inlets  and  ducts  and  turbomachinery 
blade  passages  and  cascades.  The  influence  of  these 
interactions  on  the  resulting  losses  and  downstream 
effects,  especially  with  incipient  separation,  is 
important  to  understand  and  predict  in  practice. 
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Because  of  its  convenient  analytical  form,  ac- 
curate  blended  representation  representation  of  the 
combined  Law  of  the  Wall  -  Law  of  the  Wake  behavior 
and  generality,  we  have  adopted  Walz’s  model  for 
the  incoming  turbulent  boundary  layer  upstream  of 
the  interaction.  For  the  low  Mach  number  adiabatic 
wall  conditions  appropriate  to  transonic  interac¬ 
tions,  it  may  be  satisfactorily  corrected  for  com¬ 
pressibility  effects  by  the  Eckert  Reference  Tem¬ 
perature  method  which  under  these  conditions  is,  in 
fact,  comparable  in  accuracy  to,  but  far  simpler  to 
implement  than,  the  Van  Driest  compressibility 
transformation  approach. 

Let  tt  be  Coles*  (incompressible)  Wake  Function, 

n  =  y/6  and  denote  for  convenience  R  r  .41  Re5*/ 

C<1  +  ")  <T  /T  )  1+W],T/T  *  1  +  .18  Me2  with  w  a  .  7  6 

we  we  w 

for  a  perfect  gas;  then  the  compressible  form  of 

Walz*s  composite  profile  may  be  written 


~)  n2(l-n)  -  2tt  +  2*  n2  (3-2 n) 


\1  +  R  / 


( . 215+. 655Rn)e 


subject  to  the  following  condition  linking  tt  to  Cf 
and  Red*:  c 


2  ’«  +  .213  +  en(l+R)  = 


Eqs.  (A-l)  and  (A-2)  have  the  following  desirable 
properties:  (a)  for  n  >  .10  or  so  U0/Ue  is  domin¬ 

ated  by  a  Law  of  the  Wake  behavior  which  correctly 


satisfies  both  the  outer  limit  conditions  U0/Ue  ♦  1 
and  dUo/dy  -*■  0  as  n  *  1;  (b)  on  the  other  hand,  for 
very  small  n  values,  U0  assumes  a  Law  of  the  Wall  - 
type  behavior  consisting  of  a  logarithmic  term  that 
is  exponentially  damped  out  extremely  close  to  the 
wall  into  a  linear  laminar  sublayer  profile  U/Ue- Rn 
as  n  -*■  0;  (c)  Eq.  (A-l)  may  be  differentiated 
w.r.t.  n  to  yield  an  analytical  expression  for 
dU0/dy  also,  which  proves  advantageous  in  solving 
the  middle  and  inner  deck  interaction  problems  (see 
text)  where  dMQ/dy  must  be  known  and  vanish  at  the 
boundary  layer  edge. 

The  use  of  the  incompressible  form  of  Eq.  (A-l) 
in  the  defining  integral  relations  for  and  6j* 
yields  the  following  relationship  that  links  the 
wake  parameter  to  the  resulting  incompressible 
shape  factor  =  ( 6 i*/ 0 )  : 


l  +  1.59y+  .  75n 

1  +  * 


Eqs.  (A-2)  and  (A-3)  together  with  the  defining  re¬ 
lation  for  R  enable  a  rather  general  and  convenient 
parameterization  of  the  profile  (and  hence  the  in¬ 
teraction  that  depends  on  it)  in  terms  of  three  im¬ 
portant  physical  quantities:  the  shock  strength 
(Me),  the  displacement  thickness  Reynolds  number 
Rc  6*  and  the  shape  factor  that  reflects  the  pri¬ 
or  upstream  history  of  the  incoming  boundary  layer 
including  possible  pressure  gradient  and  surface 
mass  transfer  effects.  With  these  parameters  pre¬ 
scribed,  the  aforementioned  three  equations  may  be 
solved  simultaneously  for  the  attendant  skin  fric¬ 
tion  Cf  ,  the  value  of  R  and,  if  desired,  the  tt 
value  appropriate  to  these  flow  conditions. 


Fig.  2 


(a)  Non-Separating 

Interaction  Pattern 


* 


I 


(b)  Typical  Interaction 
Zone  Properties 


NORMAL  SHOCK 


T*(y)  - 


linearized 

SUBSONIC 


DISPLACED  EDGE 


-<f. 


ROTATIONAL  INVISCID  FLO*  KITH 
FROZEN  TOTAL  SHEAR  STRESS 


f  rj 

THIN  SHEAR-DISTURBANCE  SUBLAYER 


Fig.  3  Triple-Deck  Structure  of  Interaction  Field  (Schematic) 
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